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ABSTRACT 

We present a method to confirm the planetary nature of objects in systems 
with multiple transiting exoplanet candidates. This method involves a Fourier- 
Domain analysis of the deviations in the transit times from a constant period 
that result from dynamical interactions within the system. The combination of 
observed anti-correlations in the transit times and mass constraints from dynami- 
cal stability allow us to claim the discovery of four planetary systems Kepler-25, 
Kepler-26, Kepler-27, and Kepler-28, containing eight planets and one addi- 
tional planet candidate. 

Key words: celestial mechanics; stars: individual (KIC 4349452, KIC 9757613, 
KIC 5792202, KIC 6949607); methods: data analysis; techniques: photometric 



1 INTRODUCTION 

To date, NASA's Kepler mission has announced the dis- 
covery of nearly a thousand candidate exoplanet systems 
including a significant number of systems with multiple 
transiting objects ( Borucki et al.|20l"T Steffen et al. 2010 



Latham et al. 20111. Such a large sample of exoplanet 



candidates can provide valuable insight into the nature 
of the general population of planets and planetary sys- 
tems, even if there were a sizeable fraction that are false- 
positive signals. Nevertheless, a sample of systems that 
are known with high confidence to be planetary enables 
a significant amount of additional science both for in- 
dividual objects and for the planet population generally. 
Moreover, it yields more efficient use of resources in these 
endeavors by indicating those systems where supplemen- 
tary observations will likely be most fruitful. 

Historically, planetary systems have been confirmed 
primarily by making dynamical measurements of the 
planet mass and orbital properties via radial velocity 
(RV) measurements. More recently ( Holman et al.|2010 
|Lissauer et al.|[2011a| ), this mass measurement has been 
accomplished by detailed modeling of transit timing vari- 
ations, or TTVs — the deviations from a constant period 
that result from gravitational perturbations among mul- 



tiple planets and their host star ( Agol et al. 2005 Holman 



and Murray 2005). Another contemporary development 



is the validation of planets by the process of elimination — 
excluding false positive signals by scrutinizing a wide va- 
riety of data on individual systems (e.g., BLENDER Torres 
|et al.| ( |2011[ ); |Fre"ssin et al.| ( |2011[ )). Yet detailed dynam- 
ical modeling of systems, the gathering of RV and other 
complementary data, and the synthesis of this informa- 
tion is time consuming for individual systems and is pro- 
hibitive for the entirety of the list of Kepler objects of 
interest (KOIs). The development of simpler, and more 
rapid, validation techniques is crucial for the timely ad- 
vancement of the exoplanetary science enabled by Kepler. 

In this work, we present a method to confirm the 
planetary nature of candidates in multiple transiting sys- 
tems. Complementary methods with the same goal are 



simultaneously being developed by Ford et al. (2011a I 



Fabrycky et al. (2011 1, and Lissauer et al. (2011c I. These 



methods are specifically intended to confirm planetary 
systems using a small number of assumptions and eas- 
ily implemented analysis. The method presented here 
broadly uses conservation of energy and dynamical stabil- 



ity to show that alternative (non-planetary) explanations 
of the observed photometric data are excluded with high 
confidence and that the objects must be both dynami- 
cally interacting and planetary in mass. Specifically, the 
properties of the TTV signature for the multiple candi- 
dates and the maximum allowed masses for objects in the 
systems demonstrate that exoplanet candidates in sev- 
eral Kepler systems are indeed planets. Generally speak- 
ing, it is difficult to construct astrophysical false positives 



that can mimic a multi-transiting system (Lissauer et al. 
|2011b| [Latham et al.|2011||Ragozzine and Holman|2010| >. 
Yet, the analysis presented herein does not invoke this en- 
hancement in planet likelihood, and we confirm the dis- 
covery of 8 total planets in the systems Kepler-25 (KOI 
244, KIC 4349452), Kepler-26 (KOI 250, KIC 9757613), 
Kepler-27 (KOI 841, KIC 5792202), and Kepler-28 (KOI 
870, KIC 6949607). 

This paper will proceed as follows. First, we discuss 
the essential components of the photometric data reduc- 
tion and the stellar properties (§21. The analysis method 
is discussed in [j3] The results of the analysis for the var- 
ious systems is found in S|4] We discuss the implications 
of these results in fj5] Finally, additional obervational in- 
formation and photometric diagnostics for some of the 
systems are given in lA] 



2 KEPLER DATA AND PHOTOMETRIC 
ANALYSIS 

The Kepler mission was designed to detect terrestrial- 
size planets in the habitable zone of the host star, ne- 
cessitating both a large sample size and sensitivity to 
a much larger range of orbital distances than ground- 



based surveys (Borucki et al. 20101. The instrument is 
a differential photometer with a wide (~ 100 square de- 
grees) field-of-view (FOV) that continuously and simulta- 
neously monitors the brightness of approximately 150,000 
main-sequence stars. A discussion of the characteristics 
and on-orbit performance of the instrument and space- 



craft is presented in Koch et al. (2010) 



2.1 Transit identification and data validation 

Each of the systems considered here was found using the 
Transiting Planet Search Pipeline (TPS) which identi- 
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Table 1. Excerpt of Transit Times for Kepler Transiting 
Planet Candidates 



KOI 


n 




TTV„ 








BJD-2454900 


(d) 


(d) 


244.01 




12.720355 + n 


x 86.087064 


244.01 





86.0821 


-0.0050 


0.0029 


244.01 


2 


111.5291 


0.0014 


0.0008 


244.01 


3 


124.2463 


-0.0019 


0.0011 


244.01 


4 


136.9696 


0.0011 


0.0010 


244.01 


5 


149.6898 


0.0010 


0.0008 


244.01 


6 


162.4084 


-0.0008 


0.0009 


244.01 


7 


175.1296 


0.0000 


0.0008 


244.01 


8 


187.8486 


-0.0013 


0.0011 


244.01 


9 


200.5706 


0.0003 


0.0024 


244.01 


10 


213.2911 


0.0005 


0.0005 



fies significant transit-like features, or Threshold Cross- 
ing Events (TCE) 



in the Kepler light curves (Jenkins 
et aLpOlO l. Data showing TCEs are then passed to the 
Data Validation (DV) pipeline ( |Wu et al.|20To| ). The DV 
pipeline fits a transiting planet model to the data, re- 
moves it from the light curve, and returns the result to 
TPS in an effort to find additional transit features. DV 
also completes a suite of statistical tests that are applied 
to the data after all TCEs are identified in an effort to as- 
sess the likelihood of false-positives. The results of these 
diagnostic tests are consistent with the planet interpre- 
tation and do not warn of potential pitfalls. 

After pipeline data processing and the photometry 
extraction, the time series is detrended as described in 



Ford et al. (2011b I in order to measure important stel- 
lar and planet candidate parameters. A first estimate for 
the stellar mass and radius (M 4 and R t ) is obtained by 
comparing the stellar T e g and logg values derived from 
an analysis of the stellar spectrum or from the Kepler 
Input Catalog (KIC, |Brown et aI1|2011[) to a set of CE- 



SAM (Code d'Evolution Stellaire Adaptatif et Modulaire, 
Morel (1997)) stellar evolution models computed in steps 
of 0.1 Mq for solar composition. 

With M» and R* fixed to their initial values, a transit 
fit is then computed to determine the orbital inclination, 
planetary radius, and depth of the occultation (passing 
behind the star) assuming a circular orbit. The transit 
lightcurve is modeled using the analytic expressions of 
Mandel and Agol| p002"t using non-linear limb darken- 
ing parameters derived for the Kepler bandpass (Claret 



2000 1 . The best fit model is found using a Levenberg- 



Marquardt minimization algorithm (Press et al. 19921. 
The best fit transit model is then removed from the 
lightcurve and the residuals are used to fit for the charac- 
teristics of next transiting candidate identified by TPS. 
The transit times obtained from these model fits are used 
in our analysis below and a brief excerpt is given in Ta- 
ble [l] Since these systems show visible TTVs we display 
their light curves in Figures[l][2j[3] and|4j where data are 
prepared (with TTVs corrected) as described in Holman 



et al. ( 2010 1 
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Figure 1. Top) Lightcurve and detrended lightcurve for 
Kepler-25. The small dots near the bottom of the panels in- 
dicate the locations of the planetary transits. Bottom) Phase- 
folded lightcurves (corrected for TTVs) for the two planets 
in the system. The periods of these planets are 6.24 days for 
Kepler-25b and 12.72 days for Kepler-25c. The excess scatter 
observed (particularly in the top panel) are due to anomalous, 
abrupt changes in reported flux being passed through the de- 
trending algorithm. 



2.2 Centroid analysis 

A centroid measurement, which gives the displacement of 
the flux centroid during transit compared with the flux 
centroid in data taken out of transit, is calculated for 
each candidate planet. The proper interpretation of the 
value derived for the centroid motion requires knowledge 
of the positions and brightnesses of the stars near the 
target. Motion can result in two different ways; a nearby 
star with the transit signature or a nearby star (or stars) 
that dilute the light from the target star, which has the 
transit signal. An analysis of the centroid motion can 
give important information in identifying the host of the 
transiting object and can help disentangle complicated, 
projected systems blended within the target aperture. 

One target star Kepler-25 is very bright (Kp = 10.7) 
and saturates the central pixels of its image on the pho- 
tometer. Thus, the centroid information can only identify 
displacements of the size of a pixel and is of limited use. 
The centroid of a second target, Kepler-27, does displace 
significantly from the nominal location during transit. 
The dimness of Kepler-27 (Kp = 15.9) presents a chal- 
lenge to the centroid analysis. A higher-resolution image 
of the target taken in the J band with the United King- 
dom Infrared Telescope (UKIRT) shows a faint, neighbor- 
ing star (Kp = 19.5) that is 2 arcseconds to the north- 
east of the target and that is not identified in the KIC 
(see Figure [5|. The centroid motion induced by this faint 
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Figure 2. Top) Lightcurve and detrended lightcurve for 
Kepler— 26. The small dots near the bottom of the panels in- 
dicate the locations of the planetary transits. Bottom) Phase- 
folded lightcurves (corrected for TTVs) for the two planets 
and the planet candidate in the system. The periods of these 
object are 12.28 days, 17.25 days, and 3.54 days for planets b 
and c and the candidate KOI-250.03 respectively. 



star is consistent with the interpretation that the faint 
star has constant brightness and the transit event arises 
from the target. The centroid analysis of the remainder 
of the eight planets shows no displacement beyond the 
3a radius of confusion for each target and therefore do 
not raise concern. When modeling the planet and stellar 
properties for each of these systems, additional light from 
neighboring stars has been accounted for. 



2.3 Stellar properties 

The Kepler Follow-up Observing Program (FOP) has ob- 
tained spectra of the host stars for three of the four sys- 
tems (Kepler-25, Kepler-26, and Kepler-27). The spec- 
tra were inspected for a prominent second set of spec- 
tral lines which would indicate the presence of an ad- 
ditional star. For host stars with modest SNR follow-up 
spectroscopy (roughly 20 per resolution element from the 
1.5m Tillinghast Reflector at FLWO), we fit spectra to 
a library of theoretical spectra using the tools of |Buch^| 
have et al. ( 2011 1. For host stars with high-SNR follow-up 



spectroscopy from Keck, we also report stellar parameters 
derived from a Spectroscopy Made Easy (SME) analysis 
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Figure 3. Top) Lightcurve and detrended lightcurve for 
Kepler— 27. The small dots near the bottom of the panels in- 
dicate the locations of the planetary transits. Bottom) Phase- 
folded lightcurves (corrected for TTVs) for the two planets in 
the system. The periods of these planets are 15.33 days and 
31.33 days for b and c respectively. 



(Valenti and Piskunov||1996| iFischer and Valenti||2005|) 



with mass and radius estimates generated from Yonsei- 
Yale isochrones. In cases without follow-up spectroscopy, 



we adopt the stellar parameters from the KIC (Brown 
|et al.|2011[ ). The adopted stellar parameters are given in 
Table [2] along with the estimate of contamination from 
nearby stars derived from the KIC. 



3 PLANET CONFIRMATION METHOD 

In a system with multiple planets, those planets in close 
proximity to or librating within mean-motion resonance 
(MMR) will induce the largest TTV signals on each other. 
The effects of other nonresonant or more distant planets 
will generally be smaller. In many instances, conserva- 
tion of energy implies that for a given pair of planets, 
the change in the orbital period of one planet due to a 
second planet will be met by a simultaneous change in 
the orbital period of the second planet but with oppo- 
site sign. The relative size of these changes is character- 
ized by the ratio of the masses of the objects. Thus, in 
a system where multiple planets are transiting, the pres- 
ence of anti-correlated TTVs between two objects is an 
indicator — sufficient, though not necessesary — that they 
are in the same system and are interacting. If one couples 
a stability requirement to this observation, the demon- 
stration that you are observing a planetary system as 
opposed to an astrophysical system is much simplified. 
Showing that the masses of two interacting planet candi- 
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Table 2. Key Properties of Host Stars 





KIC-ID 


Kp 


Contam^] 




log 9 


[M/H] 


v sin i 


R* 




Source^] 










(K) 


(cgs) 




km s — 1 


R® 


M 




244 (25) 


4349452 


10.734 


0.019 


6190(80) 


4.23(0.07) 


0.01(0.09) 


11.2(0.4) 


1.36(0.13) 


1.22(0.06) 


TRES 


250 (26) 


9757613 


15.473 


0.077 


4500(100) 


4.5(0.2) 


-0.21(0.08) 


1.9(0.1) 


0.59(0.03) 


0.65(0.03) 


Keck (SME) 


841 (27) 


5792202 


15.855 


0.206 


5400(60) 


5.4(0.20) 


0.41(0.04) 


0.6(5.0) 


0.59(0.15) 


0.65(0.16) 


Keck (SME) 


841 (27) 








5260(10) 


4.67(0.23) 


0.23(0.39) 


2.76(1.52) 






Keck (SPC) 


870 (28) 


6949607 


15.036 


0.090 


4590 


4.29 


0.34 




0.70 


0.75 


KIC 



a The Kepler system number is given in parenthese (e.g., 25 means Kepler-25). 

k Contamination values come from the KIC except KOI-841 which includes the effects of the non-KIC star seen in high resolution imaging. 
c Quoted uncertainties do not include systematic uncertainties due to stellar models. Values with no stated uncertainties are from KIC. 
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Figure 4. Top) Lightcurve and detrended lightcurve for 
Kepler-28. The small dots near the bottom of the panels in- 
dicate the locations of the planetary transits. Bottom) Phase- 
folded lightcurves (corrected for TTVs) for the two planets in 
the system. The periods of these planets are 5.91 days and 
8.99 days for b and c respectively. 



dates must be planetary gives enough evidence to claim 
a planet confirmation. 

When planets interact, their interactions can be 



manifest at a variety of natural frequencies (see e.g., Agol 
et al.| ( |2005| ); |Nesvorny and Morbidd!!] ( |2008| )) and the 
obital period changes at all of those frequencies may show 
anticorrelation. Nevertheless, near MMR only a few such 
frequencies will dominate the TTV signal. In the pres- 
ence of timing noise, the largest signal will be the first 
to appear with corresponding signals at other frequen- 
cies emerging over time. With this in mind, we set out 
to identify the dominant frequencies in the transit time 
residuals for all of the objects in multiply-transiting sys- 
tems. 















lit 



























Figure 5. UKIRT J-band image of Kepler— 27. A faint star 
that is not found in the KIC is located 2" from the target and 
can be seen just to the left of the target. 



After first calculating the residuals from a linear 
ephemeris (cf. Figure |7J, a sinusoidal function is fit to 
the data, using the Levenberg-Marquardt algorithm, at 
a variety of periods: 

{2wtPi\ , D f2ntPi' 

— =— + B cos - 



fi — Asm 



V p 



V p 



+ C (1) 



where A, B, and C are model parameters, P is the mean 
orbital period of the planet in question, and Pi is the 
test timescale. The fitted values for A and B are stored 
along with their measured uncertainties, a a and era, de- 
rived from the covariance matrix of the three parame- 
ter fit — the constant component (typically very small) 
is discarded, having no information at the specified fre- 
quency. This is essentially a Fourier transform of the 
(non-uniformly sampled) data. 
Next, a quantity S 

w = _ / A X A 2 B^B 2 

\ cr AlO"A2 O\B10\B2, 



(2) 



is calculated for each of the sampled periodiQ where the 
"1" and "2" subscripts correspond to the two objects. If 

For computational reasons there are small differences be- 
tween the periods (Pi) sampled in the application of equation 
[l]to each KOI as the algorithm is applied as a pipeline to data 
on all KOIs at a set of fixed periods relative to each KOI's 
orbital period. Consequently, when calculating H, the closest 
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the TTV signal from each object is anticorrelated at a 
particular period (or frequency), then S will be a large 
positive number (cf. Figure |8j. Finally, the maximum 
value that S has for a given candidate pair, H max , is the 
statistic we choose to demonstrate that the planets are 
interacting. 

Once E max has been found for a candidate pair, we 
use a Monte Carlo test to find the probability that the ob- 
served anticorrelation is caused by random fluctuations. 
To do this, we randomly shuffle the O-C residuals (and 
their errors), assigning them without replacement to the 
different transit epochs. We then refit the period and off- 
set, and apply equations |TJ and |2| to the shuffled data 
(cf. Figure [9]|. The goal here is not to maintain correla- 
tions (hence the data are not rotated, Rosary-style, when 
shuffled) but to identify the probability that data with 
the observed variance in the O-C residuals can reproduce 
or exceed the observed S max . 

Should additional correlation (or anticorrelation) in 
the O-C residuals beyond the frequency corresponding 
to H max exist, the most likely explanation is additional 
terms in the equations describing the dynamical interac- 



tions (sec, for example, Equations (A7) and (A8) in Agol 



et al. ( 2005 l) or interactions with other planets in the 
system. The generation of an anticorrelated TTV signal 
among two planets from star spots or correlated photo- 
metric noise (red noise) is virtually impossible as the dy- 
namical signal occurs over very long timescales, involves 
two different orbital periods (hence two different noise 
frequencies are required), and is only sampled on the rel- 
atively short timescales of ingress/egress and the transit 
duration. 

We run f 4 realizations of the Monte Carlo test. Any 
system where fewer than 10 random realizations have a 
S m ax greater than the nominal data is considered to be 
interacting with a false alarm probability FAP < 10 -3 . 
In the four cases confirmed in this manuscript, there were 
no examples where random data had a more significant 
anticorrelation than the original data. An extrapolation 
of the tail of the distribution of E max for each system 
allows for an estimate of the false alarm probability. Such 
estimates are generally much more significant than what 
can be obtained with our 10 4 realizations. 

Once the systems have been established as interact- 
ing via the Monte Carlo simulation, we run dynamical 
simulations to test the stability of the system. We are 
particularly interested in the maximum allowed masses 
of the objects in the system on stability grounds. Start- 



R„ 



where M v is in 



ing with nominal masses of M p 

Earth masses and R v is in Earth radii jLissauer et al. 



2011b I, we scale the masses to large values with a com- 



mon multiplier (see Fabrycky et al. (20111 for additional 
details). We integrate these systems until they become 
unstable and then reduce the mass multiplier and re-run 
the integration — iterating until the systems are long-term 



relevant periods are used when comparing objects in the same 
system. As the periods sampled initially with equation [l] are 
very dense compared to the dynamical periods of the systems 
in question, there is negligible effect on the calculation of S 
except at the shortest few periods, Pi. 
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Figure 6. Results of stability tests for the various systems. 
We adopt as the maximum allowed mass, the smallest mass 
that was observed to cause the system to be unstable. In all 
cases, the maximum allowed mass is planetary. 



stable (for ~ 10 7 years). The smallest mass that is shown 
to be unstable is chosen as the maximum allowed planet 
mass in the system. If a system is shown to be both inter- 
acting via the observed anticorrelation of the O-C resid- 
uals and has a maximum allowed mass that is planetary, 
we claim that the objects are confirmed planets. 

We note that in several cases the mass of the planets 
is very likely to be much smaller than the stated limits. 
Other considerations with bearing on mass limits, but 
which might take additional care in their application, in- 
clude mass limits based upon the observed TTV ampli- 
tudes and constraints on physically plausible densities (in 
the absence of significant and highly contrived contami- 
nation sources, some of the densities of these planets at 
the stated mass limits exceed 1000 g cm -3 ). Thus, while 
there may be some misestimation of the upper mass limits 
(e.g., caused by using a common multiplier or from uncer- 
tainties in the planet radius), there are other quantities 
that can be brought to bear on systems where stability 
does not yield sufficiently low masses. 



4 RESULTS 

4.1 General properties of the systems examined 

We first present the results of the stability test for the 
systems under study, and then show the results of the 
Monte Carlo test of the anticorrelation statistic H. Fig- 
ure [6] shows the maximum stable mass estimates for all 
of the objects in question. In all cases the masses are 
certainly substellar with only two objects being near the 
13 Jupiter mass (Mj) gray area — the rest have maxi- 
mum masses smaller than 10 Mj. As argued above, the 
actual masses of the planets in these systems are likely 
to be much smaller than the maximum masses allowed 
from stability considerations. Should one of the systems 
ultimately prove to be host to a planet/brown dwarf or 
multiple brown dwarf system, that alone would mark a 
novel discovery. 

The maximum allowed masses from this stability test 



© 0000 RAS, MNRAS 000,[2lfl4] 



Kepler planetary systems 25-28 7 



are given in Table [3] Also in that table are the values for 
the planet orbital period and transit epoch, planetary ra- 
dius, orbital distance, number of observed transits with 
their median uncertainty, and the median absolute devia- 
tion of the O-C residuals. Table[3]has similar information 
on the unconfirmed planet candidate, KOI-250.03. 



4.2 Results: Kepler-25 

The two planets in Kepler-25 have a period ratio of 2.039, 
just outside the 2:1 MMR and the O-C residuals are 
shown in Figure [7] Even if not actually resonating, one 
would expect large orbital perturbations in such a sys- 
tem. The observed anticorrelation in the TTV signal is 
indeed quite large. A plot of H is shown in Figure [8j The 
Monte Carlo analysis of the transit times shows that the 
probability of observing an anticorrelation as large as ob- 
served is much smaller than the chosen 10 -3 threshold. 
The distribution in the H max statistic for Kepler-25 is 
shown in Figure [9] Extrapolating the tail of the observed 
distribution in 5 max by eye shows that the actual false 
alarm probability is likely much smaller than 10 _3 j^] 

Dynamical stability tests, described above and in 



Fabrycky et al. (2011 1, were run on this system in order 



to establish that the objects must have planetary masses. 
The maximum allowed mass for the two planets in this 
system (shown in Figure ki| are 4.16 Mj (3.97xlO _3 M ) 
for Kepler-25c and 12.7 Mj (1.21 x 10~ 2 Af Q ) for Kepler- 
25b — meaning that a system with these masses and the 
observed orbital phases is predicted to be unstable on a 
timescale of roughly 10 7 years. 




200 300 400 

HJD -2454900 




200 300 400 

HJD -2454900 



Figure 7. Observed minus calculated transit times for 
Kepler-25. Planets b and c are the top and bottom panels 
respectively. 



4.3 Results: Kepler-26 

The two planets in Kepler-26 have a period ratio of 
1.4045, in this case slightly inside the 3:2 MMR. Here 
again one would expect large orbital perturbations. The 
Monte Carlo analysis of this system shows the probability 
of a spurious anticorrelation to be below our threshold 
(a rough extrapolation places it near 10" 5 ). The tran- 
sit times for the two planets in Kepler-26 are shown in 
Figure [10] and a plot of S is shown in Figure |11| The 
distribution in the H max statistic from the Monte Carlo 
for Kepler-26 is shown in Figure |l2"| The dynamical tests 
show that the maximum allowed mass for the two planets 
in this system are 0.380 Mj (3.63 x 10~ 4 A'/ Q ) for Kepler- 
26b and 0.375 Mj (3.58 x 1O" 4 M ) for Kepler-26c. 

An interesting feature in this system is a smaller 
peak in Figure [Tl] near a period of 90 days. Taken alone, 
this might not warrant much notice. However, a similar 
peak in S, at the same period is evident for the interaction 
between the planet Kepler-26b and the planet candidate 
KOI-250.03 and a trough in 5 is found at that period in 
the analysis of Kepler-26c and KOI-250.03 — indicating 
that they have a correlated TTV signal at that frequency 




Figure 8. Plot of S vs. period for Kepler— 25. 



instead of an anti-correlated one (see Figures 13 and 141. 
This is a hint that these three objects may be mutu- 
ally interacting such that when the middle planet's pe- 
riod increases, the periods of the two outer objects both 
decrease — and vice versa. The current analysis is not suf- 
ficient to demonstrate that KOI-250.03 is a real planet, 
but with additional data or with a more sophisticated, 
three-body analysis such a goal may well be achieved. 



2 We note that two sets of Gaussian deviates, when analyzed 
in this same manner, produce a similar structure to what is 
seen in Figure [8] between periods of 30 and 120 days, and the 
Monte Carlo generates a histogram similar to the large bump 
in Figure [9] 



4.4 Results: Kepler-27 

The period ratio of the two planets in Kepler-27 is 2.043, 
similar to Kepler-25 — slightly outside the 2:1 MMR. The 
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Table 3. Key Properties of Planets and Planet Candidates 



KOI 


Planet 




Epoclpj 

(d) 


P 

(d) 


or 


R P b ¥\ 
p U 


a 6 > c 
(AU) 


nTTpl 


(d) 


MAEpj 

(d) 


_A^p ^ yr ajj*^ j 

M Jup 


244.01 


Kepler- 


-25c 


86.0871 


12.7204 


2.89 


4.5 


0.110 


36 


0.0008 


0.0009 


4.16 


244.02 


Kepler- 


-25b 


73.5126 


6.2385 


3.58 


2.6 


0.068 


71 


0.0021 


0.0023 


12.7 


250.01 


Kepler- 


-26b 


78.8321 


12.2829 


2.82 


3.6 


0.085 


35 


0.0030 


0.0044 


0.380 


250.02 


Kepler- 


-26c 


82.8854 


17.2513 


2.12 


3.6 


0.107 


24 


0.0040 


0.0045 


0.375 


250.03 






69.2705 


3.5438 


1.98 


1.3 


0.037 


132 


0.0231 


0.0129 




841.01 


Kepler- 


-27b 


91.6726 


15.3348 


3.47 


1.0 


0.118 


30 


0.0048 


0.0069 


9.11 


841.02 


Kepler- 


-27c 


86.4274 


31.3309 


4.75 


1.9 


0.191 


15 


0.0032 


0.0029 


13.8 


870.01 


Kepler- 


-28b 


75.6227 


5.9123 


2.77 


3.6 


0.062 


76 


0.0083 


0.0057 


1.51 


870.02 


Kepler- 


-28c 


81.7277 


8.9858 


4.39 


3.4 


0.081 


52 


0.0081 


0.0102 


1.36 



a BJD-2454900 

b From Borucki et al (2011) 

c Updated to reflect stellar properties from Table [2] 
^ Number of transit times measured in QO-6 

e Median absolute deviation from linear ephemeris measured during QO-6 
f Based on assumption of dynamical stability and stellar mass from Table 2 



Table 4. Results from Monte Carlo simulations 



Figure 9. Results of the Monte Carlo test for Kepler— 25. The 
observed S max (dark gray) is significantly larger than any of 
the tested samples. 



KOI 



False Alarm Probability 



244.01 


244.02 


< 10" 3 


250.01 


250.02 


< 10~ 3 


841.01 


841.02 


< 10~ 3 


870.01 


870.02 


< 10~ 3 




Fabrycky et al. 


738.01 


738.02 


0.0007 


806.02 


806.03 


< 10~ 3 


935.01 


935.02 


0.4021 


952.01 


952.02 


< 10~ 3 




Ford et al. 




168.01 


168.02 


< 10" 3 


1102.01 


1102.02 


< 10~ 3 



Monte Carlo analysis of this system shows the probability 
of a spurious anticorrelation to be below our threshold 
(again likely to be much smaller). The transit times for 
the two planets in Kepler-27 are shown in Figure [15] and 



The distribution in the 
17| The dynamical tests 



a plot of 5 is shown in Figure [16 
H max statistic is shown in Figure 
give maximum masses for the two planets in this system 
of 9.11 Mj (8.69 x 1O" 3 M ) for Kepler-27b and 13.8 Mj 
(1.32 x 1O" 2 M ) for Kepler-27c. 



4.5 Results: Kepler-28 

In Kepler-28 the ratio of orbital periods is 1.52, out- 
side the 3:2 MMR. The Monte Carlo analysis of this sys- 
tem shows the probability of a spurious anticorrelation to 
be extremely tiny (probably of order 10 -7 ). The transit 
times for the two planets in Kepler-28 are shown in Fig- 
ure[l8]and a plot of 5 is shown in Figure [T9[ The distribu- 
tion in the H max statistic for Kepler-28 is shown in Figure 
|20| The dynamical tests give maximum masses for the two 



planets in this system of 1.51 Mj (1.44 x 1O" 3 M ) for 
Kepler-28b and 1.36 Mj (1.30 xl(r 3 M a ) for Kepler-28c. 

4.6 Other Kepler systems 

In addition to the systems presented here, similar stud- 
ies were conducted for the systems analyzed in |Fabrycky| 
et al.| ( |2011[ ); |Ford et al.| ( |2011a| >. Table [I] shows these 
results along with the results of the two other studies for 
the same planet pairs. There is broad agreement about 
the planetary nature of the sytems investigated. Never- 
theless, some differences in the results of the three meth- 
ods is an indication of their complementarity — certain 
systems, and the data we have on them, lend themselves 
better to certain types of analysis. 



5 DISCUSSION 

The Kepler mission has produced a significant number 
of new transiting exoplanet candidates, a large portion 
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HJD -2454900 

Figure 10. Observed minus calculated transit times for the 
three objects in Kepler-26, the panels are for planets b (top), 
c (middle), and KOI-250.03 (bottom). 




Figure 11. Plot of S vs. period for Kepler-26. 
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Figure 12. Results of the Monte Carlo test for Kepler-26. 




Figure 13. Correlation graph for Kepler-26b and KOI-250.03. 
The peak near 90 days is marginally significant and may indi- 
cate anticorrclated TTVs between these two objects. 



of which are found in systems with multiple transiting 
objects. One challenge presented by such a rich yield is 
confirming or validating the planetary nature of these 
candidates — especially since many of the Kepler systems 
are very dim (Kp > 15) and are very small (R < 4-R©), 
and are therefore not amenable to the traditional spec- 
troscopic methods that have been the workhorse of exo- 
planet discovery. 

In this paper (and its companion papers) we present 
a method that can fill this need. We claim that, in multi- 
transiting systems, the demonstration that two objects 
are in the same system and that their masses must be 
planetary are a sufficient criteria to confirm the objects as 
planetary. The first criterion effectively eliminates back- 
ground astrophysical false positive scenarios while the 
second shows that the objects are of planetary mass. We 
have shown that anti-correlated transit timing variations 
of sufficient significance can satisfy the first criterion. We 
developed a Fourier-based analysis and defined a statistic 
H max that can be used to calculate this significance of an 
anti-correlation observed for a pair of TTV signals. 

We use dynamical stability as our primary means 
to constrain the mass of the transiting objects. These 
tests demonstrate that the object's masses in four Ke- 
pler systems satisfy our second criterion given above. An- 
other means to constrain planetary masses comes from 
the TTV signals themselves. If the masses of the transit- 
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Figure 14. Plot of 3 vs. period for Kepler-26c and KOI- 
250.03. The dip near a period of 90 days indicates possible 
correlated (instead of anticorrelated) TTVs. Such a situation 
would arise if two planets are being simultaneously perturbed 
by a third object. The two objects whose correlation is shown 
here lie on either side of Kepler-26b, which may be the cause 
of this apparent correlation. 




200 300 400 

HJD -2454900 
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O -20 




200 300 400 

HJD-2454B00 



Figure 15. Observed minus calculated transit times for 
Kepler— 27, the top and bottom panels are for planets b and c 
respectively. 



ing objects are too large, then the associated TTV signal 
will also be too large to be consistent with the observed 
variation. 

The application of this method to Kepler data 
confirms the planetary systems Kepler-25, Kepler-26, 
Kepler-27, and Kepler-28. This method is also applied 



to the systems announced in Fabrycky et al. (20111 and 



Figure 16. Plot of H vs. period for Kepler-27. 
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Figure 17. Results of the Monte Carlo test for Kepler-27. 



the conclusions found in those studies. The fact that there 
are some differences among the three methods presented 
in this set of papers shows the importance of complemen- 
tary analysis — particularly now at the early stages of the 
application of TTV methods to large quantities of data. 

As this, and other methods continue to be refined, 
there is a significant opportunity for both the discovery of 
new planetary systems, for the elimination of astrophysi- 
cal false positive explanations for observed transit signa- 
tures, and for the identification of systems that merit ad- 
ditional, more detailed characterization. The benefits to 
exoplanetary science of transiting exoplanets, and partic- 
ularly of multi-transiting exoplanetary systems, are well 
known ( Charbonneau et al.|2007 Ragozzine and Holman 
20101. A fast and reliable method to identify exoplanet 



Ford et al. (2011a I with results broadly consistent with 



systems from large quantities of photometric data (specif- 
ically Kepler data) is an essential element for rapid ad- 
vancement in the field. Methods such as that presented 
here can inform and economize a wide variety of supple- 
mental (generally ground-based) observational efforts. 
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Figure 19. Plot of H vs. period for Kepler-28. 
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APPENDIX A: ADDITIONAL DATA 

Given the anticorrelated TTV signals in these systems 
and the maximum allowed masses of the objects, we claim 
that the objects causing the observed transit features are 
planets. Additional follow-up data from a variety of ob- 
servatories was also gathered on these systems. Here we 
show these data and their corresponding analyses for the 
different planetary systems. 




Figure Al. Ks-band adaptive optics image of Kepler— 25. The 
image is 20" on a side. 



Al High resolution imaging 

High resolution images exist for Kepler-25 and for 
Kepler-26 in the form of Adaptive Optics and Speckle 
observations. These images can be used to identify the 
presence of stars that would either cause contamination 
or be the source of a false positive system. For Kepler-25, 
Figure [AT] shows a Ks-band image that is 20" on a side 
(north is up and east is to the left) taken at the Palomar 
Observatory. 

Figure |A"2"1 is a speckle image taken at 692nm with the 
WIYN telescope. This image (and a separate image that 
is not shown) demonstrates that Kepler-25 is an isolated 
star within a radius of 0.05" for companions that are 5.5 
magnitudes fainter in R and 4.5 magnitudes fainter in 
V. Figure |A"3] is a speckle image taken at the same wave- 
length for Kepler-26. This image shows that Kepler-26 is 
a single star within a similar radius for companions that 
are fainter by 3.2 magnitudes in R and 2.8 magnitudes in 
V. 



a> 0.0 




-1.0 



1.0 0.5 0.0 -0.5 -1.0 

Distance (arcsec) 

Figure A2. Speckle image of Kepler-25 taken with a filter 
centered at a wavelength of 692nm. 



A2 Spitzer observations 

Kepler-25c, Kepler-25b, and Kepler-26c were ob- 
served during several transits with W 'arm- Spitzer /IR AC 



( Werner et aL]|2004| |Pazio et al.||2004| ) at 4.5 urn (pro- 
gram ID 60028). Each visit lasted 7 h 06 min, 8 h 30 min, 
and 5 h 45 min for Kepler-25c, Kepler-25b, and Kepler- 
26c respectively. The data were gathered in full-frame 
mode (256 x 256 pixels) with an exposure time of 12 s 
per image, yielding 1696, 2054 and 2450 images, respec- 
tively. From the images we produced a photometric time 
series using the method described in |Desert et al. (20091 
which consists of finding the centroid position of the stel- 
lar point spread function (PSF) and performing aperture 
photometry using a circular aperture on individual Basic 



Calibrated Data (BCD) images delivered by the Spitzer 
archive. 

These files are corrected for dark current, flat- 
fielding, detector non-linearity and converted into flux 
units. We converted the pixel intensities to electrons us- 
ing the information given in the detector gain and expo- 
sure time provided in the FITS headers; this facilitates 
the evaluation of the photometric errors. We adopt pho- 
tometric apertures which provide the smallest errors; the 
optimal apertures are found to be at 3.0 pixels. We find 
that the transit depths and errors vary only weakly with 
the aperture radius for all of the light-curves. Outliers in 
flux and positions greater than 5 a were rejected using 
a sliding median filter and the first half-hour of observa- 
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Figure A3. Speckle image of Kepler-26 taken with a filter 
centered at a wavelength of 692nm. 



tions, which are affected by a significant telescope jitter 
before stabilization, are alse rejected. 

We estimate the background by fitting a Gaussian 
to the central region of the histogram of counts from the 
full array. Telescope pointing drift results in fluctuations 
of the stellar centroid position, which, in combination 
with intra-pixel sensitivity variations, produces system- 
atic noise in the raw light curves. The final photometric 
measurements used in the analysis are presented in Ta- 
ble [XT] and the raw time series are presented in the top 
panels of Figures |A4] |A5j and [A7) 

We correct the light curves for instrumental effects 
and measure the transit depths and their uncertainties as 



described in Desert et al. (2011 1. The transit light curves 



model used is the 1DL transit routine OCCULTSMALL from 
|Mandel and Agol| ( [2"002| |. Only the ratio R p /R+ is allowed 
to vary while the other model parameters are set to their 
value derived from the Kepler lightcurve and the mid- 
transit times are fixed at the measured central transit 
time. We simultaneously fit the instrumental functions 
with the transit for each individual visit, and report the 
values in Table |A"T| The errors on each photometric point 
were assumed to be identical, and were set to the rms of 
the residuals of the initial best-fit obtained. 



A3 Study of possible blend scenarios 

One issue that can arise with systems like those we 
present is that, while the transiting objects are known 
to be planets, there is a small probability that they orbit 
a background or physically associated star rather than 
the Kepler target star. The centroid offsets alone can 
eliminate a significant amount of blend scenarios. How- 
ever, other types of analysis can boost confidence that 
the target star is indeed the planetary host star. 

A more detailed analysis with BLENDER folds in data 
from a number of sources and is sufficiently powerful to 
validate the planetary nature of an object on its own by 
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Figure A4. Observation of a transit of Kepler-25c with 
Spitzer. The raw and unbinned lightcurve is shown in the top 
panel. The red solid lines correspond to the best fit model, in- 
cluding the planetary transit and the time and position instru- 
mental decorrelations. The bottom panels show the corrected, 
normalized, and coadded lightcurve that is binned by 17 min- 
utes. The green curve shows the expected transit model from 
the Kepler spacecraft which agrees, within the uncertainty, 
with Spitzer. 




Figure A5. Observation of the transit of Kepler-25b with 
Spitzer. The raw and unbinned lightcurves are shown in the 
top panel. The formatting is the same as in Figure |A4] 



the elimination of possible scenarios ( |Torres et al.|2011[ ). 
However, such a study consumes a significant amount of 
human and computer resources. We conducted a small- 
scale analysis for blend scenarios for the systems pre- 
sented here, though given the expense of even these stud- 
ies we anticipate that future applications of the planet 
confirmation method described herein will not include a 
similar blend analysis. In this less exhaustive blend study, 
we quantify the likelihood that the planets we see transit 
an unresolved star in the photometric aperture whether 
by the chance alignment of a background star, or a phys- 
ical companion to the target. 

The TTV signal provides a strong constraint on the 
maximum mass ratio planets could have with the star 
they orbit. This maximum mass ratio translates to a max- 
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Table Al. Warm-Spitzer observations of KOI-244.01, KOI-244.02, and KOI-250.02 at 4.5 urn. 
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Figure A6. Observations of the transit of Kepler-26c with 
Spitzer. The formatting is the same as in Figure |A~4] 



Figure A7. Observations of the transit of Kepler-26c with 
Spitzer. The formatting is the same as in Figure \Ka\ 



imum radius allowed as a function of the spectral type of 
the blend. We use the upper envelope of the known ex- 
oplanets mass-radius diagram to estimate the maximum 
radius — from which we get the maximum dilution factor 
(or magnitude difference) per spectral type that would 
match the observed transit depth. When the maximum 
allowed magnitude difference is small, it excludes a large 
fraction of background blends as well as the reddest stel- 
lar companions. 

Stellar companions can also be constrained by the 
color of the target. We verify that simple stellar mod- 
els using T ff and the Solar isochrone are consistent with 
the r — K color (using r from the KIC and 2MASS K). 
This excludes blends that are significantly redder than 
the target (having a A(r — K) > 0.1). For Kepler-25 and 
Kepler-26, the Spitzer observations are consistent with a 
planetary transit and show no evidence for a false positive 
scenario. We exclude the scenario of a larger planet tran- 
siting a significantly redder star for both physical com- 
panions and for background stars for systems where there 
would be a 3cr discrepancy between the Spitzer and Ke- 
pler observations. 

We compute the frequency of blends using a target- 
specific (magnitude and position) Besancon model of the 
galaxy for background stars with a cut on the maximum 
magnitude difference and spectral type. For the frequency 
of physical companions, we use the rates given in Ragha- 



van et al. (20101 and draw, for our samples, the mass 



ratio between the companion and primary from the dis- 



tributions given in Duquennoy and Mayor ( 1991 1. 

For Kepler-25 and Kepler-26, we estimate that it 
is > 300 times more likely that the planets transit the 
target star than they transit an unresolved companion or 
background star. For Kepler-27 and Kepler-28, it is more 
difficult to exclude the scenario where the planets transit 
an unresolved stellar companion. The centroid shift does 
not help with this issue, leaving the maximum mass ratio 
and color as the remaining constraints. Table |A2| shows 
the results of our analysis for blend scenarios. Even for 
cases where companion stars remain a possibility, this 
does not affect the planetary nature of these systems — 
nor does it affect the system where the planets are located 
(only the star within that system). Indeed, if these plan- 
ets happen to orbit a smaller star the maximum allowed 
mass from out stability study would, in turn, be smaller 
as the dynamics cares about the ratio of the masses and 
not the actual masses; and the anticorrelated TTV signal 
already demonstrates that these are not planets orbiting 
different stars. 

This paper has been typeset from a Te^/ IMfX file pre- 
pared by the author. 
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Table A2. Probability of selected blend scenarios. 



KOI max AKp max Radius of Background Physical 
Confusion (") Blends Companion 



244 4.5 0.069 0.0000019 0.0(F] 

250 3.0 0.10 0.0000252 0.0(F 

841 4.3 0.24 0.000559 0.114 

870 5.5 0.78 0.00714 0.340 



a Strong constraints on physical companions comes from the Spitzer data on 



these targets. 
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